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Abstract 

A recently developed formalism in which Kohn-Sham calculations are combined with an "average 
pair density functional theory" is reviewed, and some new properties of the effective electron- 
electron interaction entering in this formalism are derived. A preliminary construction of a fully 
self-consitent scheme is also presented in this framework. 



I. INTRODUCTION 



Kohn-Sham (KS) Density Functional Theoryi^ (DFT) is nowadays one of the most pop- 
ular methods for electronic structure calculations both in chemistry and solid-state physics, 
thanks to its combination of low computational cost and reasonable performances. The ac- 
curacy of a KS-DFT result is limited by the approximate nature of the exchange-correlation 
energy density functional i^dn]. Typical cases in which present-day DFT fails are strongly 
correlated systems, the description of van der Waals forces, the handling of near degeneracy. 
Much effort is put nowadays in trying to improve the DFT performances via the construc- 
tion of better approximations for the KS £^ c [n] (for recent reviews, see, e.g., Refs. Q^O), 
or via alternative routes, like, e.g., the use of non-KS options.- A popular trend in the de- 
velopment of new KS -EecM is the use of the exact exchange functional E x [n] (in terms of 
the KS orbitals), and thus the search for an approximate, compatible, correlation functional 
E c [n\. 

In this work we review the basis of a theoretical framework- 1 ^ in which KS-DFT is com- 
bined with an "average pair density functional theory" (APDFT) that provides an explicit 
construction for £? c [n], transfering the work of finding an approximate functional to the 
search of an effective particle-particle interaction. A self-consitent scheme for this approach 
is presented, and some new properties of the effective interaction that enters in this combined 
formalism are derived. Very preliminary applications are discussed. 

II. DEFINITIONS 

Our target problem is finding the ground-state energy of the standard iV-electron hamilto- 
nian in the Born-Oppenheimer approximation (in Hartree atomic units, ft = m = a = e = l, 
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used throughout), 

H = T + V ee + V ne , (1) 
1 N 

T = ~Y,Vl (2) 

i=l 
1 N 1 

V « = 9 E 1 |> ( 3 ) 

AT 

V ne = 2^Wne(r,), (4) 

i=l 

where w ne is the external potential due to nuclei. Given ty, the exact ground-state wave- 
function of H, we consider two reduced quantities that fully determine, respectively, the 
expectation values 1 V^. e 1 and (\l/|Ke|^), i.e., the electronic density, 

n(r) = N ^2 J |^(r(Ti,r 2 (T 2 , ...,r N cr N )\ 2 dr 2 ...dr N , (5) 

CT1...(TJV 

and the spherically and system-averaged pair density /(r^) (APD), which is obtained by 
first considering the pair density P 2 (r!,r 2 ), 

^2(ri,r 2 ) = N(N - 1) ^ / |*(ri<7i, r 2 cx 2 , r 3 a 3 , r N a N )\ 2 dr 3 ...dr N , (6) 

a 1 ...a N 

and then by integrating it over all variables except r 12 = |rj — r 2 |, 

/(r 12 ) = l|p 2 ( ri ,r 2 )^R, (7) 

where R = |(ri + r 2 ), r 12 = r 2 — ri. The function /(r 12 ) is also known in chemistry as 
intracule density^^^^^, and, in the special case of a an electron liquid of uniform 
density n, is related to the radial pair-distribution function g{ri 2 ) by g(r) = 2f(r)/(nN). 
We thus have 

(*\Vne\*) = J n(r)v ne (r)dr (8) 
= / ^dr 12 = f -^A,r\ 2 dr l2 . (9) 

J r 12 Jo r l2 

In the following text we will also deal with modified systems in which the external po- 
tential and/or the electron-electron interaction is changed. Thus, the notation V ne and V ee 
will be used to characterize the physical system, while the modified systems will be defined 
by V = Y2iLi v { T i) an d W — \ w (\ r i ~ r il)) where the pairwise interaction w always 
depends only on [r» — r,-|. 



III. THE EXCHANGE-CORRELATION FUNCTIONAL OF KS-DFT 



In standard DFT one defines a universal functional of the one-electron density n(r) as 
resulting from a constrained search over all the antisymmetric wavefunctions \& that yield 



F[n; Vee] = mm(^\T + V ee \V), (10) 
or, more completely, clS 8b Legendre transform 16 

F[n- V ee ] = sup |min(#|T + V ee + V\V) - J n(r)«(r)dr J . (11) 

In both Eqs. (fTU|) and pip , the dependence on the electron-electron interaction has been 
explictly shown in the functional. The universality of the functional F stems exactly from 
the fact that the e-e interaction is always l/r^- The ground-state energy Eq of the system 
can then be obtained by minimizing the energy with respect to n, 

E = nun If[u; V ee \ + J n(v)v rie (v)dr^ . (12) 

A possible way to derive the Kohn-Sham equations in DFT is to define a set of hamilto- 
nians depending on a real parameter Ai^^, 

H x = T + W x + V x } (13) 

having all the same one-electron density, equal to the one of the physical system 

n x (r) = n(r) VA. (14) 

If W XphyB = V ee and W x=0 = (e.g., W x = XV ee ), one can slowly switch off the electron- 
electron interaction, while keeping the density fixed via a suitable external potential V x . 
Obviously, the APD /(r 12 ) changes with A. By the Hellmann-Feynmann theorem, 

dE x /rXl dW x dV x . rX , f j-x/ ,dw x {r 12 ) , f . ^dv x {r) , 

~d\ = W-a\ + -ax^ = J &\ 2 + J (15) 

so that by directly integrating Eq. ( fl5l) . and by combining it with Eq. ( fl2|) . one obtains 

/-^phys /■ r)iii x ( r \ 

dX J dr 12 f x (r 12 ) ^ , (16) 

where T s [n] = F[n; 0] is the kinetic energy of a noninteracting system of N spin-^ fermions 
with density n(r). The adiabatic connection in DFT thus naturally defines the Kohn-Sham 



non- interacting kinetic energy functional T a [n]. The second term in the right-hand-side 
of Eq. ffT6l) is an exact expression, in terms of the APD f x (r '12), for the Hartree and the 
exchange-correlation functional, E H [n] + E xc [n}. The one-body potential at A = is the 
familiar Kohn-Sham potential, t> A=0 (r) = t>Ks( r )- 

The traditional approach of DFT to construct approximations for £^[71] is based on the 
idea of universality. For example, the familiar local-density approximation (LDA) consists in 
transfering, in each point of space, the pair density from the uniform electron gas to obtain 
an approximation for / A (r 12 ) in Eq. (|T6|) . Our aim is to develop an alternative strategy 
in which realistic APD f x (r \ 2 ) along the DFT adiabatic connection are constructed via a 
formally exact theory that must be combined with the KS equations in a self-consitent way. 
The formal justification for this "average pair-density functional theory" (APDFT) is the 
object of the next Sec. [TV] 



IV. AVERAGE PAIR DENSITY FUNCTIONAL THEORY 

As shown by Eqs. jS]) and dHJ), the APD /(r 12 ) couples to the operator V ee in the same 
way as the electronic density n(r) couples to V ne . In order to derive an "average pair density 
functional theory" (APDFT) we thus simply repeat the steps of the previous Sec. II I II by 
switching the roles of / and n, and of V ee and V ne .- 

We thus define a system- dependent functional (i.e., a functional depending on the external 
potential V ne , and thus on the specific system) of the APD /(r 12 ) as 

G[f;V ne ] = min(*\T + V ne \V), (17) 
#-►/ 

where, again the minimum is over all antisymmetric wavefunction \l/ that yield a given 
f(r 12). We can also define the system-dependent functional G as 

G[f; V^] = sup |min(*|T + W + V ne \^>) - J f{r 12 )w{r 12 )dr 12 } . (18) 

The ground-state energy could then be obtained as 

E = min ((?[/; V ne ] + [ ^-dr 12 ) , (19) 
/eAT/ I J r 12 ) 

where Mf is the space of all A-representable APD (i.e., coming from the contraction of 
an A-particle antisymmetric wavefunction). The definition of the space Mf is evidently 
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related to the iV-representability conditions for the pair density, for which recent interesting 
progresses have been made 2 ^. In our case, however, we combine APDFT with DFT so that 
the minimization of Eq. f)19p is never directly carried on. 

In order to find the analog of the KS system for APDFT, we define an adiabatic connection 
similar to the one of Eq. (1131) in which, this time, we switch off the external potential. We 
thus introduce a set of hamiltonians depending on a real parameter £, 

H* = T + + (20) 

in which the function /(r 12 ) is kept fixed, equal to the one of the physical system, 

fHrn) = /(na) V£. (21) 

If — y ne an( j y£=° = o (e.g., = £V ne ), we are switching continuously from the phys- 

ical system, to a system of N free electrons interacting with a modified potential w^ =0 (r 12 ). 
That is, f{r\2) is kept fixed as £ varies by means of a suitable electron-electron interaction 
while the one-electron density n(r) changes with £. Again, by the Hellmann-Feynmann 
theorem, we find 

dEl , T£ ,dWt dV*,. f ... dwt{r 12 ), f € , .dv*(i 



~ = (^\^~ + ~^\^} = J f(ri2)—^dr l2 + J nZ(v)-^dv, (22) 



so that 



G[f- V ne ] = Ttif] + di J drn^v)^-, (23) 

where Tf[/] is the kinetic energy of a system of N free (zero external potential) interacting 
spin- 1 fermions having the same /(r 12 ) of the physical system. In the case of a confined 
system (atoms, molecules) the effective interaction u^ =0 (ri2) must have an attractive tail: 
the hamiltonian corresponding to £ = in Eq. (1201) describes a cluster of fermions whose 
center of mass is translationally invariant. The functional Tf [/] is the internal kinetic energy 
of this cluster. 

To fix the ideas, consider the simple case of two electrons, e.g. the He atom. When 
£ = 0, we have two fermions in a relative bound state (similar to the case of positronium, 
but with a different interaction). This relative bound state is such that the square of 
the wavefunction for the relative coordinate r\ 2 is equal to /(n.2) of the starting physical 
system. The corresponding effective interaction w^ =0 (r 12 ), obtained 7 by inversion from a 
very accurate wavefunction, 2 ! is shown in Fig. [1], for the case of the He atom. 
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FIG. 1: The electron-electron interaction at the two ends of the APDFT adiabatic connection of 
Eqs. PD11 - PI]) in the case of the He atom. 



For more than two electrons, 7f[/] is still a complicated many-body object. Moreover, 
since the corresponding i^ =0 (r 12 ) can have an attractive tail (as in the case N = 2), we may 
have an "exotic" true ground-state for the cluster, i.e., the cluster state with the same f(r \ 2 ) 
of the physical system can be an excited state. However, we have to keep in mind that our 
aim is not to solve the many-electron problem by means of APDFT alone: we want to use 
APDFT to produce realistic /(V 12 ) along the DFT adiabatic connection of Eqs. (fT3l) - (|T^I) . 
To this end, we proposed^ an approximation for the functional Tf[f] based on a geminal 
decomposition, 

T g \f\= mm V^l-V^}, (24) 

where ipi(r 12 ) are some effective geminals (orbitals for two electrons, but only depending on 
the relative distance ri 2 ), and $i some occupancy numbers to be chosen. For example, one 
can always make a "bosonic" choice, by occupying only one geminal,— equal to a/ f{ji 2 ). 
The geminals ipi(ri 2 ) then satisfy the equations 



[-V r 2 12 + w eS (r 12 )}ipi{r 12 ) = €i ipi{r 12 ) 

ZAMn 2 )\ 2 = f(n 2 ), 



with 



ST g \f) 

n 2 ' 5/(r-i2) 5f(r 12 ) 



w e s{r 12) = h 



(25) 



(26) 
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The approximation of Eq. (|2"%1) is mainly motivated by the need of having simple equations 
for f{r 12) (the one-dimensional character of Eqs. (125]) is of course very appealing). Notice 
that only in the case N = 2 we have T g [f] = T f [/], and thus the effective interaction w e fj(r 12 ) 
of Eqs. (T35]) becomes equal to u^ =0 (r 12 ). 



V. COMBINING DFT AND APDFT IN A SELF-CONSITENT WAY 

As explained in the previous sections, to compute the expectation value of the physical 
hamiltonian of Eq. ([T]) we only need /(r 12 ) and n(r) for V ee and V ne , and the non-interacting 
KS kinetic energy T s [n] plus the APD f x (r 12) along the DFT adiabatic connection (Tl~3T) for 
the expectation value of T; schematically: 

(V\H\V}= (tf|T|tf> +(y\V ee \y} + (y\V ne \V}. (27) 

T s [n]+/ A (n 2 ) /( m ) n(r) 

Our aim is to obtain n(r) and T s [n] via the KS equations, and / A (?" 12 ) via Eqs. (1251) . which can 
be generalized to any hamiltonian along the DFT adiabatic connection, by simply replacing 
the physical hamiltonian H with H x of Eqs. (|T3|) - (fl4l) in the steps of Sec. IIVI 
A self-consitent scheme for this construction reads 

(T + V KS ) $ KS = £ks$ks => n(r), T s [n] (28) 

[-V r 2 12 +^ e A ff (r 12 ;[n])]^ A (r 12 ) = e A ^ A (r 12 ) ^ 
E,^l^ A (ri 2 )| 2 = / A (r 12 ), 

E = wm{T,[n]+J o d\ J dv 12 f x (r 12 ) 9w g ^ u) + j dv n{v)v ne {v) J . (30) 

The computation starts with a trial ^Ks( r ) m the KS equations, schematically represented by 
Eq. (128]) . where $ks is the Slater determinant of KS orbitals. From the KS equations we thus 
get a first approximation for the density n(r) and the non- interacting kinetic energy T s [n]. 
Provided that we have a prescription to build an approximate w x s for a given density n(r) 
(see next Sec. lVip . we can obtain / A (r 12 ) along the DFT adiabatic connection from Eqs. (1291) . 
In general, this step is not expensive: Eqs. (129]) are unidimensional, and if the dependence 
of u> A (ri 2 ) on A is smooth, few A values (5-20) are enough to provide a good estimate of 
the coupling-constant average. The physical ground-state energy E is then evaluated via 
Eq. ( 130]) . The procedure should then be repeated by optimizing the KS potential vks, so that 
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Eq is minimum. The iV-represent ability problem of the KS exchange-correlation functional 
is clearly shifted to the iV-representability problem for f x (r 12). In view of the new conditions 
derived for the pair density,™ this seems to leave space for improvements. 

VI. PROPERTIES OF THE EFFECTIVE ELECTRON-ELECTRON INTERAC- 
TION 

So far, we have only replaced the problem of finding an approximation for i^cM with 
the problem of constructing w^ s (r 12; [n]). In order to proceed further, we study here the 
properties of Wg ff (r 12 ; [n]). 

If we want our Eqs. fl28|) - fl30|) to be fully self-consitent, we should impose that for A = 
Eqs. (|29|) yield / A=0 (r 12 ) = iks( r i2), i-e., the same APD we would obtain by inserting the 
KS Slater determinant of Eq. (I25|) in Eqs. ©-0- This corresponds, in the usual DFT 
language, to treat exchange exactly. The first property we should thus impose to w^ s (r 12 ) is 

^e ff = °(r 12 )=^ K ff S (r 1 2). (31) 

If we use only one geminal to define T g [f] in Eq. the property ( 13TT) corresponds to 

W7g i f°(r 12 ) = V 2 A/iks( r i2)/ \j iks( r i2)- For more than one geminal we need more sophisti- 
cated constructions, mathematically equivalent to those used to construct the KS potential 
^Ksfj) f° r a given spherical density n(r).— Equation (I3T]) also provides a very good starting 
point to build w^iryi): the KS system already takes into account the fermionic structure 
and part of the effect of the external potential in the physical problem. What is then left, 
that needs to be approximated, is the effect of turning on the electron-electron interaction 
without changing the one-electron density n(r), and the difference between Tf[/] and T g [f]. 

For confined systems (atoms, molecules) another property to be imposed on ^ ff (r 12 ) 
concerns the eigenvalue e^ ax corresponding to the highest occupied geminal in Eqs. ff29l) . 
In fact, the asymptotic behavior of the pair density P 2 (ri,r 2 ) of Eq. ([6]) for |ri| — > 00 (or 
I r*2 1 — » 00) is, in this case,— 

lim P 2 (ri,r 2 ) =n(r 1 )njv-i(r 2 ){ri}, (32) 

|r].|-KX3 

where riN-i(r) is one of the degenerate ground-state densities of the (N — l)-electron system 
(in the same external potential V ne ), with the choice depending parametrically upon the 
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direction = / 1 r 1 1 . A similar expansion holds for the KS pair density, obtained from the 
KS Slater determinant of Eq. ff28l) . 

lim P 2 KS ( ri ,r 2 ) = n(n)n^_ 1 (r 2 ){f 1 }, (33) 

|ri|— >oo 

where we have used the fact that, by construction, the iV-electron density is the same 
for the physical system and for the KS one (while, of course, the corresponding (N — 1)- 
electron densities are in general different). For a given attractive (atomic, molecular) external 
potential vanishing at large distances, the iV-electron density is in general more diffuse 
(decaying slower at large distances) than the (N — l)-electron density, so that the asymptotic 
behavior of the APD f(r 12) is, for large r\%, dominated by the iV-electron density decay at 
large distances. We thus see, from Eqs. (1321) and (I33l) . that the corresponding APD's, /(ri 2 ) 



and /Ks( r i2) ; wm have the same large-r 12 decay, oc e v ^^"™** ri2 , with a different prefactor 
(which can also include a polynomial function of 7*12), depending on the difference between 
n7v_i(r) and nj^ s _ 1 (r). Since the same expansion holds for any P 2 A (ri,r 2 ) along the DFT 
adiabatic connection, 

lim P 2 A ( ri ,r 2 ) =n(r 1 )n*_ 1 (r 2 ){r 1 }, (34) 
|n|->oo 

the highest eigenvalue e max in Eqs. ff29l) must be independent of A and equal to the one for 
the KS APD, 

A _ A=0 _ KS (ok) 
max max max" \ UKJ J 

In particular, if we choose only one geminal^ for the definition of T g [f], there is only one 
eigenvalue, which must be the same for every A. 

For an extended system we have scattering states in Eqs. ( 1291) . For the special case of the 
uniform electron gas, Eqs. (1291) become equivalent to an approach that was first proposed 
by Overhauser,— and further developed by other authors in the past five years.— i2L2& In 
this approach, the geminal occupancy numbers $j are the same as the ones for a Slater 
determinant: occupancy 1 for singlet states (even relative angular momentum £), and occu- 
pancy 3 for triplet states (odd relative angular momentum £), up to N(N — l)/2 geminals. 
Rather simple approximations for the effective potential w^ s (r i 2 ) gave good results^^ 1 ^ 
for the radial distribution function g(r 12 ), when compared with quantum Monte Carlo data. 
The long-range asymptotic behavior, in this case, corresponds to a phase-shift sum rule for 
the interaction w^ s (r 12 ).-2S The choice of one geminal for the uniform electron gas has been 
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explored, with remarkable success, in Ref. |30|. In this case, the formal similarity with the 
Fermi-hypernetted-chain approach 31 (FHCN) was exploited to build a good approximation 
for u> A ff (r 12 ) (which was split into jf and fj, contributions). 

Finally, the small-r 12 behavior of the effective potential u>g ff (r 12 ) is determined by the 
choice of the adiabatic connection path. For instance, if we choose u> A (r 12 ) = A/r 12 in 
Eqs. (TT3l)- (jT6{) . then the APD / A (ri 2 ) displays the electron-electron cusp / A (ri 2 — ► 0) = 
/ A (0)(l + Ar 12 + ...), which implies, in turn, that also w^ s must behave, for small r 12 , as 
^eff( r i2 — > 0) = X/r 12 + .... In particular, for A = A phys , we always have 
1/^12 + as shown, for the case of the He atom in Fig. [U If we choose a cuspless nonlinear 
path, like w A (r 12 ) = erf (Ar 12 )/r 12 , then the small r 12 behavior of u> A ff (r 12 ) is known only 
when we are approaching the physical interaction.— 

VII. PRELIMINARY APPLICATIONS 

The construction of an approximate w A ff (ri 2 ; [n]) can thus start with the decomposition 
«4(ri2; W) = w f&(ru) + w\r 12 ) + A«4(r 12 ; [n]), (36) 

where the term Aw7 A ff (r 12 ; [n]) should take care of the fact that, when the electron-electron 
interaction is turned on, the one-electron density n(r) and (for confined systems) the highest 
eigenvalue e A ax do not change. 

As a starting point, we applied the method of Eqs. (I28l) - (l30l) to the He isoelectronic 
series. In this simple (yet not trivial) 2-electron case, we have the advantage that we can 
treat Tf [/] exactly. We developed an approximation^ for Au> A ff (r 12 ; [n]) based on the one used 
for the uniform electron gas.— This approximation is designed to mimic the conservation 
of n(r) along the DFT adiabatic connection, but does not take into account the eigenvalue 
conservation. It works remarkably well when combined with a nonlinear adiabatic connection 
path w A (r 12 ) = erf (A r 12 )/r 12 that separates short- and long-range correlation effects, and 
is reported in the Appendix of Ref. |6|. Preliminary implementations of the self-consitent 
procedure of Eqs. (I28 l) -( l30l) yield ground-state energies within 1 mH with respect to full 
configuration interaction (CI) calculations in the same basis set. However, the way we carried 
out these first tests was simply based on a direct minimization of few variables parametrizing 
vxs( r )- This rather inefficient way to implement Eqs. fl28l) - fl30l) needs further improvment. 
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Besides, the Kohn-Sham potentials we obtain in this way are very unstable (like the ones of 
Ref. |33j), although the corresponding total energies and electronic densities are stable, and 
do not display the variational collapse of perturbation-theory-based approximate E c [n]M- 



VIII. PERSPECTIVES 

The generalization to many-electron systems of nonuniform density of the approximation 
built in Ref. [o for w^iru) is not straightforward. If we simply apply it to the Be atom 
case (by using only one geminal), we obtain energy errors of 300 mH. Adding the eigenvalue 
conservation can improve the results, but, of course, there is not a unique way to impose it. 
So far, we found that the final outcome strongly depends on how we impose the eigenvalue 
conservation (i.e., if we use a functional form with one parameter adjusted to keep the 
eigenvalue independent of A, the results drastically depend on the chosen functional form). It 
seems thus necessary to switch to more than one geminal, and/or to find better constructions 
for w^(r 12 ). 

In particular, it may be promising to explore the possibility to construct approxima- 
tions inspired to the FHNC,— £L2H and to try to include some of the new results on N- 
representability conditions for the pair density— Different approximations with respect to 
the one of Eq. (J21J for the functional Tf[f], and the use of Eq. (J2"3"j) also deserve further 
investigation. 

Acknowledgments 

We thank E.K.U. Gross for useful discussions and suggestions. One of the authors 
(P.G.G.) gratefully aknowledges the 30th International Worksohop on Condensed Matter 
Theory organizers for supporting her participiation to the meeting. 



1 W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

2 A.E. Mattsson, Science 298, 759 (2002). 

3 C. Fiolhais, F. Nogueira, and M. Marques (eds.), A Primer in Density Functional Theory 
(Springer- Verlag, Berlin, 2003). 



12 



4 J. P. Perdew, A. Ruzsinszky, J. Tao, V. N. Staroverov, G. E. Scuseria, and G. I. Csonka, J. 
Chem. Phys. 123, 062201 (2005). 

5 see, e.g., T. Leininger, H. Stoll, H.-J. Werner, and A. Savin, Chem. Phys. Lett. 275, 151 (1997); 
R. Pellet, A. Savin, T. Leininger, and H. Stoll, J. Chem. Phys. 116, 1250 (2002); J. Toulouse, F. 
Colonna, and A. Savin, Phys. Rev. A 70, 062505 (2004); J. Toulouse, F. Colonna, and A. Savin, 
J. Chem. Phys. 122, 014110 (2005); J. G. Angyan, I.C. Gerber, A. Savin, and J. Toulouse, Phys. 
Rev. A 72, 012510 (2005); E. Goll, H.-J. Werner and H. Stoll, Phys. Chem. Chem. Phys. 7 , 
3917 (2005); E. Goll, H.-J. Werner, H. Stoll, T. Leininger, P. Gori-Giorgi, and A. Savin, Chem. 
Phys. 329, 276 (2006). 

6 P. Gori-Giorgi and A. Savin, Phys. Rev. A 71, 032513 (2005). 

7 P. Gori-Giorgi and A. Savin, Philos. Mag. 86, 2643 (2006). 

8 A.J. Thakkar, A.N. Tripathi, and V.H. Smith, Jr., Int. J. Quantum Chem. 26, 157 (1984), and 
references therein. 

9 C.A. Coulson and A.H. Neilson, Proc. Phys. Soc. London 78, 831 (1961). 

J. Cioslowski, B.B. Stefanov, A. Tan, and C.J. Umrigar, J. Chem. Phys. 103, 6093 (1995). 

1 J. Cioslowski and G. Liu, J. Chem. Phys. 109, 8225 (1998). 

2 E. Valderrama, J.M. Ugalde, and R.J. Boyd, in Many-electron densities and reduced density 
matrices, edited by J. Cioslowski (Kluwer Academic/Plenum Publishers, New York, 2000). 

3 E.R. Davidson, Reduced Density Matrices in Quantum Chemistry (Academic Press, New York, 
1976). 

4 A.J. Coleman and V.I. Yukalov, Reduced Density Matrices: Coulson's Challenge (Springer- 
Verlag, New York, 2000). 

5 M. Levy, Proc. Natl. Acad. Sci. U.S.A. 76, 6062 (1979). 

6 E. Lieb, Int. J. Quantum Chem. 24, 243 (1983). 

7 J. Harris and R. Jones, J. Phys. F 4, 1170 (1974); D.C. Langreth and J.P. Perdew, Solid State 
Commun. 17, 1425 (1975); O. Gunnarsson and B.I. Lundqvist, Phys. Rev. B 13, 4274 (1976). 

8 W. Yang, J. Chem. Phys. 109, 10107 (1998). 

9 A. Savin, F. Colonna, and R. Pollet, Int. J. Quantum Chem. 93, 166 (2003). 

:0 P. W. Ayers and E. R. Davidson, Int. J. Quantum Chem. 106, 1487 (2006); P. W. Ayers, Phys. 

Rev. A 74, 042502 (2006); P. W. Ayers and E. R. Davidson, Adv. Chem. Phys., in press. 
:1 D.E. Freund, B.D. Huxtable, and J.D. Morgan III, Phys. Rev. A 29, 980 (1984). We used 

13 



an improved version (provided to us by C. Umrigar) of the accurate variational wavefunctions 
described in this work. See also C.J. Umrigar and X. Gonze, Phys. Rev. A 50, 3827 (1994), and 



10. 



Ref. 

22 A. Nagy, J. Chem. Phys. 125, 184104 (2006). 

23 F. Colonna and A. Savin, J. Chem. Phys. 110, 2828 (1999); Q. Zhao, R. C. Morrison, and R. 
G. Parr, Phys. Rev. A 50, 2138 (1994); R. van Leeuwen and E. J. Baerends Phys. Rev. A 49, 
2421 (1994). 

24 M. Levy, J.P Perdew, and V. Sahni, Phys. Rev. A 30, 2745 (1984); M. Ernzerhof, K. Burke, 
and J.P. Perdew, J. Chem. Phys. 105, 2798 (1996). 

25 A.W. Overhauser, Can. J. Phys. 73, 683 (1995). 

26 P. Gori-Giorgi and J.P. Perdew, Phys. Rev. B 64, 155102 (2001). 

27 B. Davoudi, M. Polini, R. Asgari, and M.P. Tosi, Phys. Rev. B 66, 075110 (2002). 

28 M. Corona, P. Gori-Giorgi, and J.P. Perdew, Phys. Rev. B 69, 045108 (2004); I. Nagy, R. Diez 
Muiho, J.I. Juaristi, and P.M. Echenique, Phys. Rev. B 69, 233105 (2004). 

29 P. Ziesche, Phys. Rev. B 67, 233102 (2003); physica status solidi (b) 242, 2051 (2005). 

30 B. Davoudi, R. Asgari, M. Polini, and M. P. Tosi Phys. Rev. B 68, 155112 (2003). 

31 see, e.g., E. Krotscheck and M. Saarela, Phys. Rep. 232, 1 (1993). 

32 P. Gori-Giorgi and A. Savin, Phys. Rev. A 73, 032506 (2006). 

33 V. N. Staroverov, G. E. Scuseria, and E. R. Davidson, J. Chem. Phys. 124, 141103 (2006). 

34 P. Mori-Snchez, Q. Wu, and W. Yang, J. Chem. Phys. 123, 062204 (2005); D. R. Rohr, O. V. 
Gritsenko, and E. J. Baerends, Chem. Phys. Lett., in press. 

35 E. Krotscheck, Phys. Lett. A 190, 201 (1994). 



14 



